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Abstract 

We extend the classical Garsia-Rodemich-Rumsey inequality to the multipa- 
rameter situation. The new inequality is applied to obtain some joint Holder 
continuity along the rectangles for fractional Brownian fields W(t,x) and for the 
solution u(t, y) of stochastic heat equation with additive white noise. 

1 Introduction 

Let the function ^ : [0, oo) —> [0, oo) be non decreasing with lim fy(u) = oo and let 

u— > oo 

the function p : [0, 1] — > [0, 1] be continuous and non decreasing with p(0) = 0. Set 
f" 1 ^) = sup<p (t)) <„ v if *(0) < u < oo 

= maxp(^)< u v if < u < p(l) 

The celebrated Garsia-Rodemich-Rumsey inequality [6] takes the following form: 
Lemma 1.1 Let f be a continuous function on [0, 1] and suppose that 



Y*^M-^>'W<B<oo. 

o Jo \ pk x — y) J 



Then for all s,t £ [0,1] we have 



r\s-t\ /4R N 

!/(*) - /(*)l < 8 / -t ]*(«)■ (!•!) 



This Garsia-Rodemich-Rumsey lemma 1.1 is very powerful in the study of the sam- 
ple path Holder continuity of a stochastic process and in other occasions. For example 



*Y. Hu is partially supported by a grant from the Simons Foundation #209206. 
Key words: Joint Holder continuity; Garsia-Rodemich-Rumsey inequality; sample path property; 
Gaussian processes, fractional Brownian fields; stochastic heat equations with additive noises. 
AMS subject classification (2010): 60G17, 26A16, 60G60, 60G15, 60H15. 



1 



if }fr(u) = \u\ p and p(u) = |u| Q+1 / p , where pa > 1, the inequality (1.1) implies the 
following Sobolev imbedding inequality 

\m - f(t)\ < C aiP \t - V> Qf 1 jT 1 '^-ff dxdy) VP . (1.2) 

The Garsia-Rodemich-Rumsey lemma has been extended to several parameter or infi- 
nite many parameters. But the parameter space are assumed to have a distance (metric 
space) and the Garsia-Rodemich-Rumsey lemma is with respect to that distance. This 
method immediately yields the following result for a fractional Brownian field W H {x) 
of Hurst parameter H = (Hi, ■ ■ ■ , Hj), then for any f3i with /3j < Hi, i = 1, • • • , d, one 
has 

d 

\W(y)-W(x)\<L^\ yi - Xi \^, (1.3) 
j=i 

where L is an integrable random variable. One can improve this result (Remark 4.5) 
by our version of multiparameter Garsia-Rodemich-Rumsey inequality. We do not seek 
for a suitable metric but rather deal directly with the multidimensional nature of the 
parameter space. 

Let us explain our motivation by considering the two parameter fractional Brownian 
field {W(x\, x 2 ), (x\, x 2 ) € [0, l] 2 } of Hurst parameter H = (Hi,H 2 ). Given two points 
x and y in R 2 , we consider the increment of W along with the rectangle determined 
by x = (xi,x 2 ) and y = (2/1,2/2): 

DW := W(yi,y 2 ) - W(x 1 ,y 2 ) ~ W(x 2 , Vl ) + W(xi,x 2 ) . (1.4) 

In [7], using a two-parameter version (1.2), the author showed that for any (3i,f3 2 with 
Pi < Hi and /3 2 < H 2 , there is an integrable random constant Lp lt p 3 such that 

\nW\<L^ 2 \yi-xi\^\y 2 -x 2 \^. (1.5) 

The above result was also obtained in [1] based on a two-parameter version of Kol- 
mogorov continuity theorem. Along the paper (in Corollary 4.4), we shall see that the 
following sharper inequality than (1.5) holds 

\UW\ < L HuH2 \yi - Xi\ Hl \y 2 - x 2 \ H ^\log(\yi - xi\\y 2 - x 2 \)\ . (1.6) 
Consequently, this estimate implies 

\W(xi,x 2 ) -W(yi,y 2 )\ < L HuH2 (\xi - yi\ H '\x 2 \ H ^\\og(\xi - yi\\x 2 \)\ 

+ \xi\ H '\x 2 - y 2 \ H ^\\og(\x 2 - y 2 \\xi\)\) 

which improves (1.3). To our best knowledge, the estimate (1.6) is the first of its kind 
in literature. We shall call such property as in (1.6) or (1.5) joint Holder continuity. It 
turns out that a large class of Gaussian fields enjoys sample path joint Holder continuity 
(Theorem 4.3.) 

Our method is first formulate and prove a multiparameter version of the classical 
Garsia-Rodemich-Rumsey inequality (1.1). The generalized inequality is then applied 
to obtain sample path joint Holder continuity for random fields. Our result generalizes 
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the results in [G], [7] and provides a different approach for sample path continuity 
problem of random fields (compare to the approach in [1], [2] and [9].) 

The paper is structured as follows. In Section 2, we shall state and prove our 
multiparameter version of the Garsia-Rodemich-Rumsey lemma. The idea is to use 
induction on the dimension of the parameter space after some observations of the 
property of operator □ defined by (1.4). Some part of the proof is similar to the 
original proof of Garsia-Rodemich-Rumsey [G] with some modification. However, we 
feel it is more appropriate to give a detailed proof. 

In Section 3, we introduce a multiparameter version of Kolmogorov continuity cri- 
teria (Theorem 3.1). To our best knowledge, a two-parameter of Theorem 3.1 first 
appeared in [1]. 

Section 4 is devoted for the study of sample path continuity for Gaussian fields. 
We give a sufficient condition for a Gaussian field to possess sample path continuity 
(Theorem 4.3). We also derive the estimate (1.6) for fractional Gaussian field. In 
Section 5, wc shall study the joint Holder continuity of solution of a stochastic heat 
equation with additive space-time white noise. 



2 Multiparameter Garsia-Rodemich-Rumsey 
inequality 

We state the following technical lemma which generalizes a crucial argument used in 
[6] in the proof of Lemma 1.1. 

Lemma 2.1 Let (Q,J-) be a measurable space and let fj, be a positive measure on 
(f2 , !F). Let g : n x [0, 1] — > M. m be a measurable function such that 

Then there exist two decreasing sequences {tk ,k = 0, 1, ■ • • } and {dk ,k = 0, 1, ■ ■ • } 
with 

tk < 4-1 = p- 1 Qp(**-i)) , * = 1,2, ■ ■ ■ (2.1) 
such that the following inequality holds 

f\q(z,t k ) — q(z,tk-i)\\ , , s 4S , . 

Proof We follow the argument in [6]. Let 

From the assumption J Q I(t)dt < B it follows that there is some to € (0, 1) such that 

I(to) < B. 

Now we can describe how to obtain the sequences dk and tk recursively for k = 1, 2, • • • . 
Given tk-i, define 



l k-l 



p- 1 Qp(f*-i)) 
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Then we choose tk < dk-x such that 

1(h) < ^- (2.3) 

and 

2^) 

It is always possible to find tk such that the inequalities (2.3) and (2.4) hold simultane- 
ously, since each of the two inequalities can be violated only on a set of tk 's of measure 
strictly less than \dk-\- Now (2.3) and (2.4) gives 

^(\ 9 M)- gM l )\\ (dz) <2£fa^ AB £ 4B 



This is (2.2). ■ 

Let x = (x±, . . . , x n ) and y = (yi, . . . , y n ) be in W 1 . We denote x' = (xi, . . . , x n _i) 
and y' = (j/i, . . . , y„_i). For each integer fc = 1, 2, ■ • • , n, we define 

Vk,yX = (Xl, . . . ,Xfc-l, J/fejXfc+l, ■ • . 

Let / be a function from M™ to M m . We define the operator 14,j/ acting on / in the 
following way: 

V kl yf(x) = f(V h ,yX) 

or Vfe,j// = / 14,j/ in short. It is straight forward to verify that 

14,1/14,?/ = 14, u 

and 

Vk,yVl,y = Vl,yVk,y 

for k =/= I. Next, we define the joint increment of a function / on an n-dimensional 
rectangle, namely 

n 

/c=i 

where / denotes the identity operator. 

Example 2.2 J/ n = 2, i/ien if is easy to see that C\ y f(x) = /(yi,j/2) — f( x i>V2) — 
f(yi,X2) + /(xi,j/2) ; which is the increment of f over the rectangle containing the two 
points x and y with all sides parallel to the axis. In particular, if f(x\ 1 X2) = x\X2, 
then \3yf(x) = (xi — yi)(x2 — 2/2)7 which is the area of the rectangle. In a more general 
case, when f has the form f(x) = IIj=i fj( x j)> then 

n 

n 

3=1 
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The following simple identity enable us to show our theorem by induction and plays 
an essential role in our approach: 



n— 1 n—l 



fe=l fe=l 
= Dyf(x',x n )-Dyf(x',y n ). (2.5) 



We are now in the position to state our general version of Lemma 1.1. 

Theorem 2.3 Let f(x) be a continuous function on [0, 1]™ and suppose that 

ID" f (x)\ \ 
* ( ' " , n K dxdy <B<oo. (2.6) 



,i]«J[o,i]" \Tlk=iPk(\xk-yk\) 
Then for all s,t £ [0, 1]™ we have 

n\si-ti\ r\s n —tn\ / |HD 

!□?/(*)! < 8' 1 / • • ■ / *-M -3 7 ) • • • ( 2 - 7 ) 



Proof We proceed by induction on n. For n = 1, it coincides with the original 
Garsia-Rodcmich-Rumsey inequality (1.1). Suppose (2.7) holds for n — 1. Let / be a 
continuous function on [0, 1]™. For any x' ,y' £ R™ -1 and any s £ [0, 1], put 



g[x,y ,s) = — ^ — 



nLift(K-^i) 

Let Q = [0, 1]™" 1 x [0, 1]™" 1 , z = (x',y')- % (2-5) we can rewrite (2.6) as 

f f U (\9M-9M\\ dzdsdt <B<oo 

Jo Jo Jn \ Pn(\s-t\) J 
Applying Lemma 2.1, we can find sequences {tk} and {dk} such that 

tk <d k -i =p- X Qpn(*fc-i)J (2.8) 



and 

p fl (z, ffc )- g (z, ffc Q|\ ^ < _45^ 



For each fceNandi'e [0, l]"" 1 , let 

f{x',t k )-f{x',t k -i) 



hk(x ) 

Again from (2.5) it follows 



Pn(|*k - ife— 1 1) 

□^/(a'.ffcJ-njrVCa/.tfc-i) 



g{x',y',t k ) - g(x',y',t k -i) 
PnQtk ~ tk-l\) 
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Thus, the inequality (2.9) becomes 



/ / 

./[O,!]™- 1 J[0, 



* \ — dx'dy' < -5 — . 



Now, by our induction hypothesis, for every k > 1, a, b € [0, 1]™ , a = (oi, . . . , 
and 6 = (b u . . . ,6 n _i), 

/.|oi-bi| /.|a„_i-6„_i| / ^no 

pr^fcWi < s"- 1 / ... / 



dpi(ui) ■ ■ ■ dpn-x{u n -i) . 

Denoting ^4 = [0, |oi— 6i|]x- • -x [0, |a n _i— 6 n _i|] and dp(ui, ■ ■ ■ , it„-i) = dpi(ui) ■ ■ ■ dp n _i(it„_i) 
the above inequality can be rewritten as 



On the other hand, by (2.8), we have 



4 n B 



v 2 d 2 

u n-l u k-l 



dp(u ir -- ,u n _i)jD n (tfc_i-ifc). 
(2.10) 



Pre(*fc— 1 — *fc) < Pn(tk-l) 

= 2p(d fc _i)<4[p(rf fc _i)-p(dfc)] 



Combining this inequality with (2.10) yields 

oo 



fc=l 

oo 



< 8"- 1 >>[p(d fc _ 1 )-p(d fc )] / 4 " 2 g )dp( Ul ,--.,un-i) 



^f^w^)-^)] /" W — 

fe=l ^ V "! 

00 ^ rdk-i r / 4™5 \ 

< V4 / / f" 1 — j ) dp(ui, ■ ■ ■ ,u n -i)dp n (u n ) 

b _ 1 Jrffc \ u l ' ' ' U n-l U nJ 



k=l 

m-iA I / ,t,-i ' 4 -° 



< 8 n_i 4 / / *"M — 2 dp( Ul , • • • , u„_i)dp„(u n ). 

With f(x' ,1 — x n ) replaced /(x', x n ) we can obtain the same bound for 

inrV&^-nrVtM)!. 

Hence, for every a, 6 G [0, l]™ -1 , 

4 n B 



- — 

A \ u l ' ' ' ' 



dp(ui,--- ,u n -i)dp n (u n ). 

(2.11) 



To obtain (2.7) for general s,t in [0, 1]™, we set 

f{t', r) = f(t', s n + T(t n - s n )) for T G [0, 1] 
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and 

p n (u) =p n (u\s n - t n \). 

Upon restricting the range of the integration in (2.6) and carrying out a change of 
variables we get 

, * B^Ml ) dxdy < 

[0,1]™ J [0,1]" \U k =l Pk{\xk - VkDPnQXn - Vn\) ) \Sn ~ t n \ 

Thus, by (2.11), we deduce 

\u n j{t)\ = p^Kt'^-n^Kt'M 

< 8" f 1 [ lsi ~ tl{ ... f^- 1 ^- 11 ^ ( i n B 



y u\ ■ ■ ■ U„\s n — t n \ 2 



dp(ui, • • • , u n -i)dp n {u n \s n -t n \). 
Another change of variables yields (2.7). ■ 

3 Sample path Holder continuity of random fields 

In this section, given a random field W, we study sample path continuity property. 
The first application of Theorem 2.3 is the following criteria for joint continuity of 
sample paths which is similar to Kolmogorov continuity theorem, which we shall call 
joint Kolmogorov continuity theorem. 

Theorem 3.1 Let W be a random field on R™. Suppose there exist positive constants 
ct-, Pk (1 < k < n) and K such that for every x, y in [0, 1]™, 



[\U n y W{x)\ a ] <K \{\x k -y k \ 



•^k yK\ 
k=l 



Then, for every e = (ei, . . . , e„) with < e^a < (3k (1 < k < n), there exist a random 
variable rj with E?/" < K , such that the following inequality holds almost surely 



\n?w(s)\<c v (u>)l[\t k -s 



k=l 



for all s,t in [0, 1]™, where C is a constant defined by 

fc=i 



C = 8"4 n/Q Y[ 1 1 



Proof Let *(w) = \u\ a , p k (u) = \up k where j k € (§, 2±^), 1 < k < n. A direct 
application of Theorem 2.3 gives that for all s, t in [0, 1]™ 



7fc - i I Mo,i] 2 ™ \Ak=i Ffc " Vk\ aik 



k=l 



(3.1) 
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m-ff *( ' n; T wl , )**• 

Mo,l] 2 " Vllfc=lPfe(^fe - Vk)J 

From our assumption and Fubini-Tonelli's theorem, 

= [[ J^ W{XT dxdy 

< K If TT K - 2/ fe | 1+ ^- tt ^^ < oo. 

Hence, the event fT = : < 00} has probability one. Therefore for each w in 

O*, the inequality (3.1) gives 

|n?W(t,o;)| < 8 n TT — — ^4 1 (4"-B(w))° 

for every s,t in [0, 1]™. For each fc, the power 7/. — — can be made arbitrarily close to 
This completes the proof with 77 = i? 1 /". ■ 

Remark 3.2 TTie result obtained by RaVchenko [7] was the inequality (3.1) in the case 
n = 2. 



4 Sample path continuity of Gaussian fields 

We now focus on sample path continuity of Gaussian random fields. In case of Gaussian 
processes (n = 1), one of the first sufficient and necessary conditions for sample path 
continuity was given by Fernique [3] (see also [5]). Namely, let p(u) be an increasing 
positive function such that 

E\W(x)-W(y)\ 2 <p 2 (\x-y\) 

for any pair (x,y) in [0,1] 2 . Then Fernique [3] showed that a sufficient condition for 
almost sure continuity of the process (W(x), < x < 1) is 

f 1 P( u ) j 

/ au < 00. 

Jo V log « 

In the original paper of Garsia-Rodemich-Rumsey [6], the authors also observed that 
the above condition is equivalent to the condition (by integration by part) 




log —dp(u) < 00. 



Later, it was shown that the above condition is also necessary [4, 8]. In case of Gaussian 
fields, recent progress on modulus of continuity of Gaussian random fields has been 
reported in [9, 2]. 
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Let W be a centered Gaussian random field with covariance function 

E[W(x)W(y)] = Q(x,y). (4.1) 

We will always assume that Q is a continuous function of x and y. For any fixed 
x, y, the random variable 0™W(x) is also Gaussian with mean zero. In the following 
proposition, we compute its variance. 

Proposition 4.1 Let W be a centered Gaussian random field with covariance function 
given by (4-1). Then 

E[\n;W(x)\ 2 ]=n 2 £ y) Q(x,x). (4.2) 
Furthermore, if the covariance function Q has the following product form 



Q(*,y) = Y[Qk(xk,y k ) (4.3) 

k=l 

then (4.2) is simplified as 

n 

E [\U n y W{x)\ 2 ] = J] [Q k {x k ,x k ) - Q k {x k ,y k ) - Q k (y k ,x k ) + Q k (y k ,y k )} . (4.4) 
fc=i 

Proof We calculate the variance directly as follows 

e [|n; i iy(a;)| 2 ] = E[n$w(x)-n$w(x)] 



= E 



nft tV) w(x)w(x) 

U^ v) E[W{x)W{x)\ 

In 

(y,y) 



The identity (4.2) follows. To prove (4.4), we notice that the pair of operators (/ 
Vk,(y,y))(I ~ Vn+k,(y,y)) transforms the k-th factor of Q in (4.3) to 

Qk(x k ,x k ) ~ Q k {x k ,y k ) - Q k (y k , x k ) + Q k (y k ,y k )- 

Since the operators I — Vk t VtV ), (1 < k < 2n) are commutative, we can write 



Ufl y) Q{x,x) = f[(I-V kt{yty) )(I-V n+k>{y , y) )Q(x,x) 

k=l 
n 

= \ [Qk{xk,x k ) - Qk(xk,yk) - Qk(yk,x k ) + Q k (yk,yk)} ■ 



k=l 



Hence, the identity (4.4) follows. ■ 

Definition 4.2 Let f be a continuous function on K™. We call a set of non-negative 
even functions {p\, . . . ,p n } joint modulus of continuity of f if 

(i) For each 1 < k < n, p k (0) = 0, and p k is non- decreasing and continuous. 
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(ii) For every pair (s,t) in 3R 2n ? the following inequality holds 

n 

pr/wi < l[Pk(\tk-s k \). 

*!=i 

In view of Theorem 2.3 and Theorem 3.1, the continuity of sample paths is governed 
by the joint modulus of continuity of Hn y sQ(x, x). Such modulus of continuity always 
exists. For instance, we can define a joint modulus of continuity for y ^Q(x,x) as 
follows. We set 



Pi(u) = sup 

x,y£[0,l] n :\xi— yi\<u 



Given pi, . . . ,Pfe-i, define 



sup 



nf« y) Q(x,x) 



x,ye[0,l] n --\x k -y k \<u H i= lpj(\Xj - Vj\) 



in which we have adopted the convention 0/0 = 0. It follows immediately that pk's are 
non-decreasing and continuous. Furthermore, we have Pfc(0) = and 



n$ ty) Q(x,x)< l[pl(\x k -y k \). 



(4.5) 



k=l 



Namely, {pi,pi,P2,P2, ■ ■ ■ ,Pn,Pn} is a modulus of continuity for Un y y jQ(x, x). We also 
call {pi, • • • ,p„} a modulus of continuity for y \Q(x, x). 

In the following theorem, we give a sufficient condition for almost sure continuity 
of a Gaussian random field. 

Theorem 4.3 Let W be a centered Gaussian random field with covariance function 
given by (4-1), and pk (1 < k < n) be a modulus of continuity for y ^Q (x , x) , namely 
the inequality (4-5) is satisfied. Suppose that 



E 

k=i 



log - dp k {u) < oo. 



(4.6) 



Then, with probability one W has continuous sample path. Furthermore, there is a 
random variable B(oj), with finite expectation, such that 



\niW{x)\ < 2-8" 



\xi-yi\ 



l0 g „2...„2 



dpx(u\)---dp n {u n ). (4.7) 



Proof We set ^(x) = e^ 2 / 4 and 
B(w) =H 



cxp 



\Oy-W(x)\'< 



m^plQxk-ykl) 



dxdy. 
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Lemma 2.3 immediately gives us (4.7) for u such that B(oj) is finite. To see (4.7) 
indeed holds for almost every uj, it is sufficient to show that B has finite expectation. 
We notice that the random variable 



N 



D"W(x) 



Uk=lPk(\xk ~ Vk\ 



is Gaussian, has mean zero and variance less than or equal to one. Thus, an application 
of Stirling's formula gives 



-exp 



EN 



2/v 



fc=0 



4 fc fc! 



Hence 



EB 



< 1 



.cxp 



(2fc)! 
=J 8* (fc!) : 

1 OO 

2 ^ 

fc=i 

iV 2 



-(EiY") 



15 
14' 



[0,1] 



dxdy < 



15 
14 



and the proof is complete. ■ 

As an application of the above theorem, we obtain joint continuity for sample paths 
of fractional Brownian field, as mentioned in (1.6). 

Corollary 4.4 Let W H be a fractional Brownian field on R™ with Hurst parameter 
H = (H\, . . . , H n ). Then, there exist an integrable random constant C n ,Hi_.--- ,n n such 
that for every x,y in [0, 1/2]™ 



p;h"*(i)| < Ohh„...,«, 



log n \ xk ~ yk \ 



\k=l 



\ \xk~Vk\ 



H, 



(4.8) 



k=l 



Proof The covariance function of a fractional Brownian field is given by 

n 

E [W H (x)W H (y)] = J] R k (x k ,y k ), 

where 



fc=i 



R k (s,t) = ^[\s\ 2H " + \t\ 2H "-\s-t\ 2H "] , VMel. 
By Proposition 4.1, we obtain the second moment for U l yW H {x) 

n 

E\a n y W H (x)\ 2 = Y[\x k ~y k \ 2H ". 

k=l 

This means that Pi(u) = u Hi , i = 1, 2, • • • , n are the modulus of continuity of ,Q(x, x). 
From Theorem 4.3 it follows that 



\xi-yi\ 



log- 



\a n y w(x)\ < c n , Hl ,.. . Hn 

Jo Jo V ui ■ ■ ■ u r , 

Now an elementary computation yields the inequality (4.8). 



\x n —y n \ 



1 



.H,-l 



■ 1 dui ■ ■ ■ du n . 
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Remark 4.5 In [2], it is shown that for fractional Brownian field W H on W 1 , there 
exist a random variable A\ = A\(uj) > of finite moments of any order and an event 
fij of probability 1 such that for any u £ O J , 



sup 



\W h {s,lu) - W H (t,u>)\ 



e[oTi]« E" =1 \sj - tj\ H * Vlog (3 + | Sj -^l" 1 ) 



Our method leads to a slightly stronger estimate. Indeed, for every s,t in [0,1/2]", 
with s and (0,0, . . . , 0, t n ), the estimate (4-8) gives the increment along an edge of the 
n-dimensional rectangle [si,ii] x • ■ ■ x [s n ,t n ] 



\W H {s u -- - ,s n )-W H (s lr - - ,s n - u t n ) 

/n-l 



k 



Kk=l 



n-l 

log n sfc 

fe=l 



1/2 



|s n - t n \ Hn | log \s n - t n 



,1/2 



Similarly, we can obtain analogue estimates along any edge of the n-dimensional rect- 
angle [si,t\] x • • ■ x [s n ,t n ]. The increment along the diagonal is majorized by the total 
increments along all the edges connecting s and t. Hence, this argument yields the 
following estimate 



\W H (s)-W H (t)\<A 1 (u;)J2 II s ? 

k=l \frk 



log n 



1/2 



Isfc-tfcl^lloglsfe-tfeH 1 / 2 . 

(4.9) 



5 Stochastic heat equations with additive space time 
white noise 

In this section let us consider the following one dimensional stochastic differential equa- 
tion 



at 



±Au + W 



«(0,iO = 



< t < T, ye 



y e 



(5.1) 



where Ai 
Pt(y) = - 



'hi 



■u, W is space time standard Brownian sheet, and W 



lo-yW- Let 



'2nt 



Then the (mild) solution of the above equation is given by 



u(t,y) 



Pt-r(y - z)W(dr, dz) . 



where the above integral is the usual (Ito) stochastic integral (however, the integrand 
is simple. It is a deterministic function). The solution u(t,y) is a Gaussian random 
field. It is known that u(t,y) is Holder continuous of exponent |— for time parameter 
and ^— for space parameter. Namely, for any a < 1/4 and any f3 < 1/2, there is a 
random constant C a $ such that 



\u(t, y) - u(s, x)\ < C a ,p (\t - s\ a + \x- yf) 



(5.2) 
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We are interested in the joint Holder continuity of the solution u(t, y). We need the 
following simple technical lemma. 

Lemma 5.1 Let a, b, 5 be some positive numbers, where a < b, and let I, J be the 
integrations 



b 



1 



1 = -= 1 - dr 



J = I —= (l-e-^)dr. 



Then for every a € [0, 1/2], 

2(Vb-Va) (l - e - ^) <I< 2{Vb-Ja) (l - e _ £) 

and 

Proof On the interval a < r < b, we have 1 — e - ^ < 1 — eT^ < I — e~^ . The 
estimate for / is then a straightforward consequence. To estimate J, we first use 
integration by part to obtain 

— / s 2 \ r=a f a ( d 5 2 \ 

J = 2Vf(l-e~-J l—e-^\2Vrdr 

2Ja~(l-e-£\+5 2 f e~^r^/ 2 dr. 



By a change of variable x = -7= , we see that 



(2 \ r°° 
l-e-fej + 2V2S / e~ x 'dx 



If -4= > 1, since lim^oo -J* — e - 2 . — 0. J is majorized by 



J < cJa~ ( 1 - e 2c 



If -4= < 1, the integration e 11 c?x is bounded by y/n/2, thus J is majorized by 



J < 2y/a(l- ) +cS 



Therefore, for any < a < 1/2, employing the elementary inequality 1 — e x < c a x a , 
we obtain 

J < c a b 2a o}l 2 - a 

and the lemma follows. ■ 

Theorem 5.2 For every a in [0, 1/4], there is an integrable random constant C a such 
that for all (s,x),(t,y) in [0,1/2] 2 



\u{t 1 y)-u{t 1 x)-u{s,y)+u{ S ,x)\<C a \t~s\ 1 / i - a \x~y\ 2a ^\\og{\t-s\\x-y\)\. (5.3) 
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Proof u(t,y) is a mean zero Gaussian field. The covariance of u(t,y) and u(s,x) is 
given by 



E[w(s, x)u(t, y)} 



X[o,s] ( r )X[o,t] (r)ps~r(x ~ z)p t - r (y - z)drdz 
f(s,x)f(t,y)drdz , 



where f(s,x) = X[o,s](r)p s -r(x - z). 

We calculate the second moment of □? x \u(t, y) as follows 



E 



v) 



®0 2 (SiX) u(t,y)n 2 (StX) u(t,y) 

En ( S!;!; , s ,x) u (^y) w (^j/) 



Oi x . x) V{u(t,y)u(t,y) 



= □ 



f(t,y) 2 drdz 



n( s ,x,,,x) [f(t,y)f(t,y)]drdz 

n 2 {S!X) f(t,y)] \n 2 StX) f(t,y) 



drdz 



where 



□(.,«)/(«. i/) 



= [f( s ,x)-f(t,x)-f(s,y) + f(t,y)Y 

= f( S ,x) 2 + f(t 7 x) 2 + f( S7 y) 2 + f(t,y) 2 

-2f(s, x)f(t, x) - 2f(s, x)f(s, y) + 2f(s, x)f(t, y) 
+2f(t, x)f(s, y) - 2f(t, x)f(t, y) - 2f(s, y)f(t, y). 
Taking the integration with respect to z and using the following identity 



p a (z - x)p b (z - y)dz = p a +b(x - y) 



we obtain 



E 



v) 

[2X[0, S ](^)P2s-2r(0) + 2X[0, t ]{r)P2t-2r(0)] dr 

+ / [-2x[o. S At]Ps+i-2r(0) - 2x[o. s ]P2 S -2r(x - y) + 2x[o, sA t]Ps+t-2r (x ~ y)] dr 
J& 

+ J [2X[0, S At]Ps+t-2r(x -y)- 2X[Q,t]P2t-2r{x ~ V) ~ 2X[0,sAt]Ps+t-2r (0)] dr 
' [P2s-2r(0) - P2s-2r{x - y)] dr + 2 / [p 2 t-2r(0) - P2t-2r{x ~ y)] dr 

o Jo 



[p s +t-2r(0) - Ps+t-2r(x - y)} dr. 
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By change of variables u = 2s — 2r, v = 2t — 2r and w = s + t — 2r in the above 
corresponding integrals respectively and noticing that s + t — 2(s A t) = \t — s\, we get 



E 



n 2 



2s i-lt 



n ( s ,x) u (^y) = / [P«(°) -Pu(x-y)]du+ / [p„(0) -jj„(x- y)]di> 
J Jo Jo 

re+i 

-2/ [p w (0) - p w (x - y)]dw 

J\a-t\ 

2(aVt) fS+t /-|s-t|\ 

- / +2 / [p r (0)-jJr(a:-l/)]dr. 

s+t J2(sAt) JO J 

By Lemma 5.1, we see that 

f2(«Vt) /-s+t \ 

[p r {0) -p r {x - j/)] dr 

< —= (l - e-fe&J (\/2s + a/27 - 2>/s+l) < 0. 



s+t J2{sf\t 



and 



Jo 

for every a in [0, 1/2]. Thus 



3 



a <Ca|a:-y| a >-*| 1/2 - a . 



An application of Theorem 4.3 immediately gives the desired result. ■ 

Remark 5.3 Using the method in Remark 4-5, the above result implies the following 
estimate 

\u(s,x) - u(t,y)\ < C (\s - t|V4^ log _i_T + \ x _ y \i/2^i og _J_Z^ (5.4) 
which is sharper than (5.2). 
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